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Abstract 

We consider the Newtonian limit of modified theories of gravity that include 
inverse powers of the curvature in the action in order to explain the cosmic ac- 
celeration. It has been shown that the simplest models of this kind are in conflict 
with observations at the solar system level. In this letter we point out that when 
one adds to the action inverse powers of curvature invariants that do not vanish 
for the Schwarzschild geometry one generically recovers an acceptable Newto- 
nian limit at small distances. Gravity is however modified at large distances. We 
compute the first correction to the Newtonian potential in a quite general class 
of models. The characteristic distance entering in these modifications is of the 
order of lOpc for the Sun and of the order of 10 2 kpc for a galaxy. 
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It has been suggested that the origin of the accelerated expansion of the universe 
might be a modification of the dynamics of gravity at low curvatures, as opposed to 
a cosmological constant, or to a cosmologically evolving scalar field (quintessence). In 
jl!2j a correction to the action of the form 5C = —Mpfj, 4 /R was proposed. M p is the 
Planck mass, R is the curvature scalar and \i is a parameter with dimensions of mass. 
It was shown that this model has cosmologically interesting solutions that can explain 
the acceleration of the universe without recourse to dark energy if /i ~ Hq. However, 
it was soon realized that this model is in conflict with solar system experiments [3] and 
presents instabilities when matter is present [3]. Several modifications were suggested 
to revive this explanation of the cosmic acceleration, like adding positive powers of 
the curvature or taking into account quantum effects to cancel the instabilities |5|6j . 
trading this correction for a singular ln(R) term in the action [7j or working in the 
Palatini formalism, that yields second order field equations instead of fourth order ones 
jH]. However one can generally expect problems in reproducing a correct Newtonian 
limit for models involving inverse powers of the curvature scalar (or Ricci tensor) only 

pro] . 

In this letter we point out that if we add to the action inverse powers of curvature 
invariants that do not vanish for the Schwarzschild solution (so they have to involve 
the Weyl or Riemann tensors) one can make an expansion of the exact solution over 
the Schwarzschild geometry and it is easy to see that at small distances the corrections 
are suppressed. This can be seen already at the level of the Newtonian action obtained 
by a weak field expansion of the relativistic one. Lets consider for definiteness a term 
in the action of the form 5C = M 2 /i 4n+2 / Q n where Q = R flua \R tJ-ucrX and n is some 
positive integer number (this kind of corrections have been considered in ^T] in the 
context of cosmological evolution). In the Schwarzschild geometry, outside a spherical 
mass distribution, this invariant reads Q = 48 ^<f^ where M is the mass of the object 
sourcing the gravitational field and r is the distance to the center of the object. We 
can then expect that the modifications induced by adding some inverse powers of this 
term to the action will be suppressed at small distances for the same reason its effect 
on cosmology becomes important only at late times, when the mean curvature of the 
universe is R ~ fi 2 . To get an order of magnitude of the expected characteristic distance 
controlling the departure from Newtonian gravity we can evaluate the Newtonian action 
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in this geometry (so we consider the weak field limit r 3> GM) 

S = J d 4 x [(V0) 2 + 8rrG(f)p\ , (1) 

where ~ —(1 + goo)/2 is the Newtonian potential, p is the energy density and = 
_QM_ f or a static point-like mass. We can now consider the shift in the action induced 
by the correction 

5S = f d A x p An+2 /Q n ~ f d A x /i 4ri+2 ^" NO , (2) 

when evaluated at the previous solution. We see that this correction will be negligible 
at small distances and will become important roughly when ^ 4 " +2 48 n(^2f) 2 " ~ (^0) 2 = 
^P~. This implies a characteristic scale for the induced modification of Newtonian 
gravity given by 

„6n+4 _ _ 

pi- 
One might think that gravity is well measured at large distances, so this kind of mod- 
ifications would be in conflict with observations, but due to the extremely small value 
of \x necessary to explain the cosmic acceleration (/i ~ H ) the corresponding distance 
at which this correction becomes important is astronomically large. For instance for 
a star like the Sun this distance is of order ~ 10 — 100 parsec, much bigger than the 
solar system, while for a typical galaxy like the milky way r c is of order ~ 10 2 — 10 3 
kilo-parsec, with lower values corresponding to lower values of n (we are taking n > 1). 
This last value is not far from the size of the galaxy itself, so this suggests that the 
modifications of gravity induced by these terms will be negligible at the level of solar 
system experiments but can have a significant impact in the dynamics of larger objects 
like a galaxy. 

To confirm this intuition we will solve at first order the modified Einstein equations, 
for a class of models, in an expansion over the Schwarzschild geometry. We will obtain 
the first correction to the gravitational potential and we will check that the scale 
controlling the departure from Newtonian gravity at large distances is indeed given by 
the formula above. So we start with the action 

S = J d A x^-g-^ [R + f(R, P, Q)] , (4) 

with 

P EE (5) 



and Q denned before as the contraction of two Riemann tensors. We are using the 
notation of [TT] . and we will also take the same form for f(R, P, Q), 

u 4n+2 

With this choice for f(R,P,Q), corrections to standard cosmology will only occur at 
the present epoch if we take /i ~ H and it has been shown that for appropriate values 
of the parameters there are cosmologically interesting accelerating solutions that are 
attractors at late times [TT] . 

Our strategy will be then to take the ansatz 

ds 2 = -[1- — + eA(r)}dt 2 + [1 - — + eB{r)]~ l dr 2 + r 2 dVL 2 2 , (7) 

and treat e as a small expansion parameter. We can write the equations of motion as 

G^ + fi 4n+2 H, u = 0, (8) 

where G^ v is the usual Einstein tensor and /i 4n+2 if Ml , is the extra term generated by 
f(R, P, Q). We expand now both tensors in e: 

G^u = G^2 + tG^2 + e G^2 + . . . , 

H,., = H^ + eH^J+e 2 Hj 2 J + .... (9) 

Since the Schwarzschild solution solves the ordinary Einstein equations, we have G^} = 
0. Treating /i 4n+2 as an order e parameter, at first order in our expansion the equations 
for A and B become (from now on we set e = 1) 

G (i) = _^ +2jf (o) _ (1Q) 



For the tt component of this equation we find: 

(2GM - r)(B + rB') ///\2+4n r e»+i 1 
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while the rr component reads: 

-2GMA + rB + r(-2GM + r)A' _ / fi\ 2 + 4 « r 6n+2 (1 + n) [An - ^-{2 + I2n) i 
r(r - 2GM) 2 ~ \2J (GM) 2n+l (3c) n (r - 2GM) 



Notice that H$ — > for r — > (and so do A and B) and we approach the Schwarzschild 
solution in this limit. Also it is clear that we can not take the limit c — > 0, since 
H$ diverges in this case and our expansion breaks down. In fact our equations are 
completely independent of a and b at the order we are working, since both R and P 
vanish in our zero-th order background. We can solve exactly the previous equations 
yielding 
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(13) 



A(r) 



f) 



4n+2 



2(1 + n)r 



6rt+3 



GM 



(6n + 1) - 4n 



(14) 



(6n + 3)(3c)"(GM) 2 ™+ 1 
This solution automatically solves the other equations because of the Bianchi identities: 
A M G^ = A fl H' MU = 0. For r ^> GM we can approximate 
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B(r) = -a 



6n + 3 GM 
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(15) 
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where a _ ( 6n+3 ) 2 4n+2 (•;,., 
potential then reads 



- and r c is defined in eq. (jSJ • The correction to the Newtonian 
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So we find indeed that the correction to Newtonian gravity becomes sizable at dis- 
tances r ~ r c . By construction, our expansion breaks down at this same distance. 
We can, however, trust our approximate solution when r r c . We can check the 
self-consistency of our expansion by inserting our approximate solution in the exact 
equation, for which we get that G^ + ix 4n+2 H^ ~ (GM/r 3 )(r/r c ) 12n+8 while the terms 
we have solved for are of order Gjl = —fi 4n+2 H^ ~ (GM/r 3 )(r/r c ) 6n+4: . In fact 
we can consider our approximate solution as the first order of the Taylor expansion in 
powers of /i 4n+2 oc l/r^ n+4 , that is generated by solving the equations of motion, order 
by order. For r ^> GM this expansion takes the form 

6n+4 
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and the validity of the expansion is controlled by the parameter r/r c . 

The full equations of motion for the metric are now fourth order, so Birkhoff's 
Theorem does not hold and the solution for a spherically symmetric mass distribution 
will not be unique in general. However only first derivatives of the metric appear in 
the first order of our expansion, so the solution is unique at this order 1 . Also, we 
have considered an expansion over an asymptotically flat geometry. One might be 
worried because flat space is not a solution in our theory, and the equations actually 
become divergent for this geometry (as r — > oo in our zero-th order background). But 
considering the spacetime asymptotically flat or not (we could consider our expansion 
in the background of one of the cosmological solutions obtained in [TT], for instance) 
would only introduce appreciable corrections at distances much bigger than r c , where 
the curvature induced by the spherical mass distribution has fallen down to the cosmo- 
logical value Q ~ Hq. And, as we said before, our expansion is only valid for distances 
r <C r c . We should mention here that our solution differs from the solution obtained 
in for the linearized version of the same theory over a maximally symmetric space. 
The reason is that the linearized expansion of the theory considered in [13] breaks down 
for small radius in the spacetime of a spherically symmetric mass, where the curvature 
is substantially larger than that of the cosmic background. 

It is remarkable that with the values of the parameters necessary to explain the 

cosmic acceleration (yU ~ H ) the correction (fTo^) is negligible for solar system scales 

(within the solar system the relative corrections with respect to the Newton potential 

are smaller than ~ iQ-io-38n^ k ecome s important at galactic lengths. So despite 

having a negligible effect on the solar system this kind of modifications would have 

implications for the Dark Matter problem, changing the required abundances of Dark 

Matter and would in fact influence the measurement of most cosmological parameters. 

It is tempting to speculate on the construction of a MOND-like modification of gravity 

[T2] along these lines that could explain the galactic rotation curves without recurring 

to Dark Matter. The model we have analyzed goes in that direction, in the sense 

that the gravitational attraction becomes stronger at long distances, but the size of 

the correction appears to be too small. Notice however that given the perturbative 

method we have used to solve the equations we can not say that much about the 

^n this case the two boundary conditions determine the mass and an overall normalization of the 
time coordinate. 
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domain r ~ r c where the extra contribution to the gravitational force is of the order of 
the Newtonian one. We leave the issue of determining the nature of the corrections at 
these scales for future work. Also, we have not studied other aspects of the model like 
the possible appearance of ghosts |13|14j and instabilities |4j. These issues are more 
model dependent 2 , and in this letter we wanted to focus on the general mechanism for 
reproducing a correct Newtonian limit at small distances in this class of models. But 
given that we are introducing a correction in the action that gets smaller as r — > 
for a spherical mass, one might also expect that any pathological behaviour associated 
with this term will only manifest itself at large distances. 
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